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Abstract 

The evaporation of the black holes during the very early universe is studied. 
Starting from black hole filled universe, the distributions of particle species are 
calculated and showed, that they differ remarkably from the corresponding equi- 
librium distribution. This may have great impact to the physics of the very 
early universe. Also the evolution of the universe during the evaporation has 
been studied. 
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It is well known that in the very early universe the particles could not been equi- 
librated due to a too short equilibration time [|l], §, §]. Indeed, the only possibility to 
have equilibrium distributions at times t ~ 10 -34 s corresponding equilibration tem- 
perature T eq ~ 10 15 K is to assume, that there is some process producing particles 
readily to equilibrium distribution. (It should be noted that in this case it is not really 
question of thermodynamical equilibrium because there are no interactions maintain- 
ing it.) Sometimes it is assumed that e.g. black hole evaporation could be such a 
process || . When the particles are originated from small black holes near Planck scale 
which evaporate due to Hawking process^], the radiation at any given time is thermal 
with temperature equal to Hawking temperature. 

In the present paper we calculate the distribution of particles originated from black 
holes. We suppose that the black holes are generated at Planck time through quantum 
gravitation. After Planck time they begin to decay emitting particles. We show that 
the resulting particle distributions are not the equilibrium one, but differs crucially 
from it. This difference may play remarkable role e.g. in inflation ||[], because the 
effective, ensemble corrected potential || differs now from the one calculated with an 
equilibrium thermal bath [0]. 

In the present study, for simplicity, we neglect the curvature effects, arising at 
Planck scale due to ambiguous concept of vacuum state ||. By assuming isotropy, 
the comoving observers have same vacuum and hence, whenever the average energy 
per particle is clearly smaller than Planck mass, the analysis should be applicable. 
Consistently we omit also all direct couplings of quantum fields to curvature and treat 
the gravitational field classically. 

Let the evolution of the universe be given by the scale parameter R(t) and suppose 
that at Planck time no other matter exists in the universe than black holes. The radius 
of a (non-rotating) black hole with mass itlbh is given by tbh = 2ms#/Mpj, where 
Mpi is the Planck mass. The black hole surface area is now given by 



For simplicity, suppose that all black holes have same mass mi at initial time and 

their number density t{ is n(ti). Because there is no other type of matter than black 

holes, the universe is at matter dominated, i.e. R(t) oc t 2 / 3 . The ratio of the 

density of the universe to the critical density at that time is now given by Q(ti) = 

87rmin(tj)/(3Mp 1 iy(tj) 2 ), where the Hubble parameter is H{ti) = ~t^~ . Thus we can 

parametrize the initial density with density ratio 

Because a black hole radiates at any time thermal radiation with the temperature 

given by the Hawking temperature Tbh = M^ x / {RimiBu) through its surface, the 

1 Thc existence of small black holes at early times would have also many other interesting conse- 
quences, see II. 



loss of energy in time unit is ABHh*Tg H 7r 2 /30, where h* is an effective number of 
particles lighter than black hole mass. So, as well known, the dynamics of a black hole 
evaporation is governed by the equation 

mBH ~ -15 x 837r»mV (2) 

where we have omitted the absorption of radiation back to the black holes. This is well 
motived whenever the radiation density £ ra d is small, i.e. Erad^BH "C |t71bh|, which is 
the case. The solution of the equation (g) reads 

( *-*A 1/3 

m B H{t) —mill , (3) 



where the black hole life-time r is given by 



5 x 8 3 7r 2 mf 



(4) 



We express the initial black hole mass in terms of the inverse of r 

M P1 / K V /3 m 

where 5 = l/(Mpir) has to be a small number <C 1, i.e. the black hole life-time is 
much longer than Planck timeP|. Using the expression for f2j and (|5|) we can give the 
black hole number density 



"( ( ')-lfUJ (6) 

We are now ready to calculate the distribution of the particles emitted by the 
black holes. We further simplify the situation by stating that all emitting particles 
are essentially massless, so that the effective particle number h* remains constant 
during the evaporation^. Strictly speaking, this is not true, but gives a reasonable 
approximation of the situation, where no specific model with its particle contents is 
given. Also, we suppose that the black holes decay before the thermalization starts. 
This means that no thermalizating processes are effective during the evaporation and 
hence not taken into account in the calculation. Whether this is a good assumption 
or not is studied later in this paper. 

2 Indced, one has to assume that U + t 3> tpi, because otherwise quantum gravitational effects 
would certainly not be omitted. 

3 The massive case is as well calculable as the massless one, but leads to introducing of numerous 
new parameters, which is avoided here. 



During a time [t, t + dt] (t > U) the energy emitted by a black hole and distributed 
to h* particle species (degrees of freedom) is dE = —rriBHit) dt. Because at a given time 
the radiation emitted obeys thermal spectrum with temperature T BH (Pbh = 1/Tbh), 
the energy per degree of freedom within momentum between p and p + dp is given by 

+ dp - dE{p) = — — — — — , n , ± ,, n — TTyT, (7) 
K J d s kuj(k)f ± {k, pBH\t)) 

where the fraction k±dE/h* represents the (total) energy per degree of freedom (k± is 1 
for bosons and | for fermions appearing due to different statistics). The other fraction 
expresses the part of the energy within the given momentum range and the functions 
f^ 1 are the usual Fermi-Dirac (+) and Bose-Einstein (-) distribution functions. Thus, 
the particle number density (at the given momentum range) is simply given by 

dE{p + dp) - dE{p) = k ± dE fHPsHp) , , 

W(p) K Jd3kLu(k)f±(f3 B Hk)- [ ' 

Taking into account that the black hole density reads n{t) = n(ti) (-R(tj) / R(t)) 3 , one 
obtains the distribution of a degree of freedom: 

d 3 p dE(p + dp) -dEjp) 

dfmMj^ = nit) ^ . (9) 

Combining @ and @ we write 

Af k±m BH (t)dt (2TT) 3 f ± ((3 BH p) 

dfm = ~ n(t) K WT(M' ( 0) 

Eq. (0) represents the distribution generated during a time-interval [t, t + dt]. 
The overall distribution thus if we want to end up to at a given time t is obtained by 
summing all contribution from ti to t, i.e. all times t', t{ < t' < t. These contribution 
are, however, red-shifted by a factor R(t')/R(t), hence the final particle distribution 
at time t is given by 



The formula above, Eq. (|Tl|), describes in principle completely the particle distribu- 
tions created by black hole evaporation. It can easily be generalized to the case, where 
the (initial) black holes are not equally massive, but have some distribution. In this 
case one just have to sum over all contributions emerging due to different black hole 
masses. 

To be able to calculate the distributions in practise, one has to solve the evolution 
of the scale parameter R(t). For that purpose we have to determine first the time t E q 



when the universe enters the radiation dominated era after the beginning of the black 
hole evaporation. We have to calculate when the energy density of the black holes, Sbhi 
equals the energy density of radiation, £ ra d, they have emitted: EBHitEQ) = SraditECj)- 
Now, the (average) energy density of the black holes at a given time t is given by 
^BH(t) = mBH(t)n(ti)(R(ti) I R(t)) 3 . During a matter dominated time R(t) oc t 2 ^ 3 and 
thus 2 

£ BH {t) = m B H{t)n{U) ( -j J , U<t< t EQ . (12) 



On the other hand the radiation energy density at t < tEQ is given by 

/d 3 k f d 3 k r d 3 k 

—-^f^^Uiik) =g f J -^—^f f (k,t)uj(k) + g b J j—^f b (k,t)cv(k), 

(13) 

where we have summed over all relevant bosonic (<7&) and fermionic (gf) degrees of 
freedom^. Eq. fll3|) can simplyfied to 

^=-{^'( ! '>® 4 *>(l!)' (i4) 

where —mBuit 1 ) (Jr^)) corresponds the energy radiated at t' by black holes and red- 
shifted until time t, whereas n{tj) (^tI) represents the number density of decaying 
black holes at the time when the radiation was produced. Evidently, Eq. (|13"D can be 
obtained directly without using any knowledge about particle distributions. Indeed, 
the connection between Eqs. fli"3| ) and QT4] ) requires the identification = g*, where 
g* is the usual effective number of degrees of freedom g* = gb + ^9f- in the rest of the 
paper we apply this identification. 

It should be noted here, that the correct evolution of the cosmic scale parameter 
at matter dominated era is not really as simple as R(t) oc t 2 ^ 3 but more complicated. 
This is due to the fact that the correct energy density of matter, i.e. black holes at 
that time reads 

£ BH (t) = m(t)n(U)(R(U)/R(t)) 3 , (15) 

which is not as simple as just dilution of the matter. Applying this equation to the 
Einstein equation, one obtains 

2/3 

(16) 
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1 (m(t)\ 



7/2 



which, however, is in quite good agreement with the usual matter dominated form. 
Moreover, the calculation of the scale parameter at the radiation dominated era leads 

4 Hcre we remind that, for sake of simplicity, we have assumed that all relevant particles are 
essentially massless, thus oj(k) = k. 



to a complicated integro-differential equation. Therefore we choose here to use the 
conventional forms of the evolution parameter R, which anyway gives good impression 
of the physics involved. It is now simple matter to take the expression for Sbh and 
£ ra d from Eqs. fli~2|) and (|14D (and relation R(t) oc t 2//3 ) and equate them. One obtains 
a rather simply equation 



m BH eq)1?eq = ~ [ dtt 2/3 m BH (t). 



(17) 



It is noteworthy, that eq. (|T7|) does not depend on the initial black hole number density 
n(ti) but only on initial time and black hole mass. According to the approximated 
analysis, the radiation dominance begins quite late compared to the black hole life- 
time. When ti/r varies from to 0.5, the combination (tsQ — U)/t decreases from 
0.928 to 0.908. In any case over 90% of the black hole life-time belongs to the matter 
dominated era. This is a consequence of the fact that the soft radiation produced at 
the early phase of the black hole life-time further red-shifts during the expansion of 
the universe. Using the expression ([IB]) for cosmic scale parameter, one obtains (tEQ ~ 
U)/t = 0.924 which is in good agreement with the approximated result. Therefore, it 
is reasonable to use the simpler equations. 

The scale factor is now simply given by equation 



R(t) 



1/2 



9{t-t 
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t 



EQ, 



+ 9{t EQ - 1) 




R(ti 



tEQ_\ 
U J 



2/3 



where 9 is the usual Heaviside-function. Eq. ( |T8"D can be inserted to formulas like (|I~T|) 
and (0). In particular we are now interested in the energy density of the universe at 
the end of the black hole evaporation era: the energy density at that time determines 
the later evolution of the universe. We obtain the energy density at a given time t 

3 

n(ti 



8{t) 



x 




R(t f ) 

dt'rfiBH(t') 



+ 9 (tEQ - t' 




(19) 



where ti <t < tf. For times t > tf the energy density is only red-shifted as usual. We 
define a pseudo-temperature T to describe the energy density and particle distribution 
by setting simply 



7T~ 

30 



g*f{t) 4 = S(t). 



(20) 



Now, T is by no means a real temperature, simply because there exists no thermal 
equilibrium. However, it describes the energy density in familiar way and, moreover, 



when the thermal equilibrium is finally maintained, T coincides with the true tem- 
perature. Therefore, at the time tf it is enlightening to compare the actual particle 
distribution to the thermal one at pseudo-temperature T. In the Figure |l]a we have 
displayed the pseudo-temperature T at tf as a function of U/t (ti = tpi). 




Figure 1: (a) The pseudo-temperature T at t = tf as a function of the ratio 
ij/r. {ti = tpi, f2j = 1.) (b) Evolution of the pseudo-temperature T as a function 
of (t — ti)/r. The upper curve corresponds U/t = 0.35 whereas the lower one 

u/t = o.i. (a = i.) 

An interesting feature emerges in the behaviour of T (i.e. energy density). For 
ti/r < z c = 0.334 the pseudo-temperature has a minimum during the black hole evap- 
oration era. In the Fig. |l|b we have plotted the evolution of two pseudo-temperatures, 
one with ti/r = 0.1 < z c (lower curve) and another with ti/r = 0.35 > z c (upper 
curve). We have scaled the time so, that T is presented as function of (t — U)/r. Note, 
that the large ti/r values correspond small black holes with evaporation times compa- 
rable to ti. Thus, if ti — tpi, the approximations are not very good and. Practically 
in all relevant cases the pseudo-temperature has a minimum. 

In this stage we have to make a notion about the homogeneity of the matter and 
radiation in the universe. All calculations done as far, and all calculations to be done, 
assume that the universe is homogeneous at large scales. We can consider the matter 
to be homogeneously distributed if there are several black holes within a single horizon 
volume. On the other hand the radiation can be viewed to be roughly homogeneous, 
if the distance between the radiation sources, i.e. black holes is shorter than the 
length propagated by radiation. These two conditions are essentially equivalent leading 



to n(t) > t 3 . In particular, this requirement at transition time t B Q means, that 
5 < 0.14(f2j) 3 / 2 . So, if fij = 1 we obtain 5 < 0.14, which in terms of black hole mass 
reads m BH > 0.356(/i,/160) 1 / 3 . 

One can also ask, when the last emitted radiation (at t — tf) is homogenized. Spa- 
tially homogeneous distribution is reached at time t^ (assuming the radiation domi- 
nation after tf) if n(t z)* -1 ' 3 R{th) I R{t /) < t^ — tf. Performing a numerical calculation 
it shows up (with Oj = 1) , that (th — tf)/r varies from to 3 as 5 increases from to 
0.5. Thus the homogenization time is of the same order of magnitude as tf itself or 
even shorter. 

Keeping in mind all above, we are now finally ready to calculate the distributions 
themselves. Applying Eqs. (0), (^), (Pf) and formulae for black hole density and 
temperature to the Eq. (JTTJ) , one obtains the particle distributions. The distributions 
can be written as 

-Wu" n °»^ [W)) f±(! wf M (21) 

In the Fig. 0a we have displayed the distribution at tf for bosons multiplied by p 2 as 
a function of the momentum p together with a reference distribution. The reference 
distribution is the equilibrium distribution with T = T(tf). All are given in Plank units 
with Qi = 1. The shape of distribution given in Fig. (^)a is generic and it is noteworthy 
that the black hole originated distribution has more power in large momenta. On the 
other hand, the IR end of the distribution, p ~ 0, there is less power. Indeed, the IR 
behaviour of reads fb ~ k/p like for equilibrium (where /~ ~ T/p), but 

K = -^mk dtmBH{t) \m) < (/) (22) 

for all relevant 5: the ratio T(tf)/n is for all 5 < 0.5 larger than ~ 5.6. This may 
be important, because the IR effects are essential in many physical considerations, in 
particular when some kind of effective potential and/or interactions are used. 

The fermionic distribution function is presented in Fig. 0b. As for bosonic case we 
have presented for comparison the equilibrium distribution at T(tf), too. The shape 
of the distribution in Fig. ^jp is generic, and again there is more power at large p in ff 
than in equilibrium distribution. In general, the distributions arising from black hole 
evaporation are flatter and lower than corresponding equilibrium distributions. This 
can be understood, because at late times t ~ tf the black hole radiates particles with 
large momenta. The particle density, h<6(/), is thus always lower than at equilibrium 
with same energy density, nffjy The ratio < 0.36 for all 5 < 0.5 and is 

proportional to when 6^0. 

Finally the question, how small the parameter 5 can be, i.e. how large the black 
holes can be. If one supposes that the black holes evaporate before the thermalization 



(a) (b) 

Figure 2: (a) Boson distribution function p 2 fb{p\tf) as a function of momen- 
tum p (solid line). For reference, also the equilibrium distribution at the tem- 
perature T(tf) has been presented (dashed line). 5 = 0.1. (b) Corresponding 
fermion distribution function p 2 ff(p;tf) (solid line) with equilibrium reference 
curve at T(tf) (dashed line). 5 = 0.1. 

takes place due to GUT processes strong interactions or electroweak interactions 
[0], one has to require that tf -C tth — where Y is the (average) thermalization 
rate of the particles involved. If T ~ aT, where inspired by GUT's a ~ 3 x 10~ 2 , we 
conclude that S/T ^> a 2 or 5 ^> 3 x 10~ 8 . This estimate has many uncertainties e.g. 
because as dimensional parameter the pseudo-temperature T has been used. Never- 
theless, this approximation gives a general idea about the time-scales involved. On 
the other hand, one should address some interest to the question, when the masses are 
negligible. Concerning the tree-level masses, general remarks are impossible without 
specifying a particular model. However, in the case of the ensemble masses ("thermal 
masses") some remarks can be done. The ensemble mass m e ~ gT, where g is a 
generic coupling, should be small compared to the black hole temperature T BH , i.e. 
m e Tbh- This can be rewritten as 5 1 / 3 /T ^> gn^-^Y^Mpl , where Tbh = Tbh(U) 
has been used because it corresponds the lowest black hole temperature giving a con- 
servative limit. Now the left hand side of this equation is always larger than 7.9 (at 
S = 0.5). Thus (for h* = 160) the ensemble mass can be ignored, if the coupling 
g « 1.7! 

In the present, paper a general study of the particle distributions inspired by the 
black hole evaporation is presented. It is found that the distribution produced differ 
clearly from the equilibrium distributions. The distribution has much more power 



at large momenta and the boson distribution has slower increase when p — * 0+. The 
conditions of validity of the analysis have been studied and found that under wide range 
of parameters calculation can be performed with good reliability. Also the evolution 
of the universe has been studied during the black hole evaporation era. In particular, 
the time when it makes transition from matter domination to radiation domination 
and the generic bahaviour of radiation density are calculated. These considerations 
may have great impact to the calculations of the very early universe. However, more 
detailed studies are definitely needed, when any particular model is considered. 

References 

[1] J. Ellis and G. Steigman, Phys. Lett. B89 (1980) 186. 

[2] K. Enqvist and J. Sirkka, Phys. Lett. B314 (1993) 298. 

[3] K. Enqvist, P. Elmfors, I. Vilja, Phys. Lett. B326 (1994) 37. 

[4] S. Hawking, Nature 248 (1974) 30; Commun. Math. Phys. 43 (1975) 199. 

[5] For review, see: E. Kolb and M. Turner, The Early Universe, Addison- Wesley, 
1990. 

[6] B.J. Carr, in Observatonal and theoretical aspects of relativistic astrophysics and 
cosmology, eds. J.L. Sanz and L.J. Giocoechea, World Scientific, Singapore, 1985; 
Ap. J. 201 (1975) 1; ibid. 206 (1976) 8. 

[7] A. Riotto and I. Vilja, Mod. Phys. Lett. All (1996) 1151. 

[8] See, e.g.: N.D. Birrel and P.C.W. Davies, Quantum Fields in Curved space, Cam- 
bridge University Press, 1982. 



